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ABSTRACT 


Recent developments in the field of flywheel powered electrical 
energy storage systems has prompted the need for a better understanding 
of the varied design and analytical criteria that must be considered in the 
selection of a flywheel. Techniques have been developed for presenting 
an analytical and graphical means for selecting an optimum flywheel sys- 
tem design, based on system requirements, geometric constraints and 
weight limitations. The techniques for creating an analytical solution are 
formulated from energy and structural principals. The resulting flywheel 
design relates stress and strain pattern distribution, operating speeds, 
geometry, and specific energy levels. The design techniques incorporate 
the lowest stressed flywheel for any particular application and achieve 
the highest specific energy per unit flywheel weight possible. Stress and 
strain contour mapping and sectional profile plotting reflect the results 
of the structural behavior manifested under rotating conditions. This 
approach toward flyw^heel design is applicable to any metal flywheel, and 
permits the selection of the flywheel design to be based solely on the 
criteria of the system requirements that must be met, those that must 
be optimized, and those system parameters that may be permitted to vary. 


INTRODUCTION 


This paper describes a procedure for designing an optimum flywheel 
shape based on the constraints of geometry, speed and stress so as to 
maximize energy density. The design procedure described relies on the 
application of linear elastic structural mechanics and the laws of con- 
servation of energy and momentum. Little work has been reported in 
maximizing the energy density of solid disc flywheels. Much work however 
has gone into the design of turbine blades and discs, and electric gene- 
rators and motors, which are perhaps the closest entity to the energy 
storage flywheel. The basic structural laws under which flywheels, tur- 
bine blades, generators and motors behave are the same but their func- 
tions, based on different design objectives, are different. 

The energy density of a flyv/heel is represented by the simple 
relationship: 



( 1 ) 
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where 


E = energy density 

Kg = flywheel shape factor (dimensionless) 

a r material working stress 

P r iTidterial density 


Flywheel shape factors for several geometries are shown in Table 1. 
For disc flywheels the shape factor can approach 1.00. The disc shaped 
flywheel that produces this high a shape factor has constant stresses 
throughout. This is attributed to the fact that each unit volume of material 
is stressed equally to a predetermined working stress level and therefore 
produces the largest amount of energy possible. The flywheel shape that 
produces this constant stress is exponential in profile. 

Equation (2) expresses the summation of forces in a flywheel 
(Reference 2). 


d(XYa^) 

dX 


. Yo +- W^X^Y 
t g 
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For uniform strength the tangential and radial stresses must be 
equal and of constant value throughout. 

Therefore , 


a* = a_ = a = constant 
t r 

Equation (2) can be restated as 


d(XY) 

dx 


0 2 2 

Y + £-W X Y = 0 
g<? 


and by integrating, 

, Y 

In = - 
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Applying the boundary conditions, at X = Xj^, Y = Y 



Substituting Equation (5) into Equation (4), 
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where 



radius 

thickness at radius X 
hub thickness 
tip thickness 


= tip radius 

Equation (6) gives the normalized thickness as a function of the normalized 
radius and the hub-to-tip thickness ratio. This equation represents the 
profile configuration for the constant stress flywheel geometry. According 
to Equation (6), the disc, even though with infinitely decreasing thickness, 
is prolonged to infinity. But practically, the disc is limited by a cylin- 
drical boundary or radius X„ at which it has a thickness Y„. Although 
the theoretical flywheel of infinite diameter would have a shape factor 
of 1.00, the practical flywheel of finite radius Xp^ would have a shape 
factor less than 1. 00, In order to improve the shape factor of this ex- 
ponentially shaped finite diameter flywheel, the author has chosen to take 
some of the material that theoretically existed between the finite diameter 
and infinity and place it near the rim of the flywheel producing a constant 
thickness section running from a point on the surface to the rim. The 
utilization of a flat tip as the means for improving the flywheel shape 
factor, in lieu of an exponentially flaired tip, was chosen for two reasons: 
number one, the shape factor difference between a flaired and a flat tip 
was found to be insignificant, and secondly, the mamufacturing and machin- 
ing operations are considerably simplfied by having a flat tip rather than 
a flaired tip. The question of how much constant thickness material 
should be added to the flywheel must now be determined. Since it is the 
objective to improve the shape factor as much as possible, it is necessary 
to solve Eqxiation (1) in terms of the shape factor, 1^, for the various 
stresses and energy densities associated with each flatted tip flywheel 
that is generated for each hub-to-tip thickness ratio used in Equation (6). 
An analytical evaluation must first be performed to evaluate the resulting 
stresses and energy densities for each flywheel geometry. 


One such analytical method developed at Lockheed Missiles Ic Space 
Company utilizes a computer program based on two dimensional stresses 
and strains developed in rotating machinery. These relationships were 
then expanded (Reference 3) and culminated in the computer program. 

A typical set of results are shown in Table 2. Tangential stresses, 
ra^al stresses and flywheel thicknesses are presented for various radii 
starting at the rim and extending to the hub. Once the stresses are deter- 
mined the program then utilizes Equation (1) to calculate the shape factor 
from the calculated kinetic energy and maximum flywheel stress. A maxi- 
mum shape factor is then obtained for each hub-to-tip thickness ratio by itera- 
tively evaluating different flatted tips that begin at different percent radii. 

Figure 1 represents the results of the relationship between the fly- 
wheel shape factor and the hub-to-tip ratio. The point at which the optimum 
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flat begins has an optimum value which varies with the hub-to-tip thickness 
ratio. By applying the appropriate flat tip as the means for optimizing the 
flywheel shape factor for a particular hub-to-tip ratio, it can be recognized 
that for hub-to-tip ratios less than 1.00, the appropriate flat tip would ex- 
tend from the tip to the hub, and the flywheel shape would be that of a flat 
unpierced disc having a shape factor of 0.606. 


FLYWHEEL PARAMETRIC TRANSFORMATIONS 


By increasing the hub-to-tip ratio, the shape factor is improved and 
becomes 1 . 00 at a hub-to-tip ratio of infinity. There are instances when 
a high hub-to-tip ratio may not be practical, such as when geometric con- 
straints limit the axial length of the hub. It is therefore desirable to 
determine the effects of changing parameters on the rest of the system. 

A set of parametric relationships for relating flywheel diameters, speeds, 
weights, kinetic energy levels, operating stresses and thicknesses permits 
an easy determination of the effects on each parameter when one or two 
are changed. Equations (7), (8), and (9) express these relationships, and 
can be used for the flywheel shapes generated by Equation (6) having fixed 
hub-to-tip ratios. 



THK 


Twr 


o ■ ('^) (^) ( 5 ) 


WT„ /THK„WD^\2 


where a = working stress 

D s flywheel diameter 
W = flywheel speed 
K£ s kinetic energy 
THK = flywheel thickness 
WT = flywheel weight 


(7) 

( 8 ) 
(9) 
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A specific example is used to demonstrate the application of 
these equations for a flywheel having a kinetic energy of 12 kilowatt-hours. 
This is shown in Figure 2 for several flywheels from three to four feet 
in diameter. The stagger of the data points is caused by the quantizing 
error in selection of either 75 or 80 percent flat value for the optimum 
hub-to-tip ratio. Since the optimum flat lies between these two values a 
smooth curve in actuality joins the optimum points. 

The curves represent a family of varying diameters for a particular 
kinetic energy level. It is reasonable to assume that the kinetic energy 
requirements are already known for a desired application. The require- 
ments for the operating speed, or the maximum flyv'heel stress will 
further restrict the number of available flywheel design selections. If we 
permit the flywheel speed to vary, we can superimpose a family of fly- 
wheel speeds on Figure 2 indicative of a specific kinetic energy and operat- 
ing stress level. This was done for an operating stress level of 100 ksi. 

By utilizing a series of plots similar to Figure 2, reflecting various operat- 
ing stress levels, a more complete selection of flywheel geometries is 
possible. These curves represent a means for selecting an optimum 
geometry flywheel based on the kinetic energy requirements, volumetric 
limitations, and desired flywheel life (reflected through operating stress 
level). Once a selection is made, a stress profile may be performed 
using the two dimensional stress program of Reference 3, which produces 
results similar to those of Table 2. If upon examination of the results of 
this initial computer run, the flywheel selected is found to be satisfactory, 
a much more rigorous, three dimensional, stress -strain examination can 
be performed using a finite element computer program. A process of con- 
tour mapping of the stresses and strains developed under rotating con- 
ditions for each and every point within the flywheel can then be made. 
Figures 3 through 6 show contours for radial, tangential, axial, and axial 
shear stresses for a quarter section view of a 12 kilowatt-hour flywheel. 
Figures 7 and 8 show the radial, tangential, axial and axial shear stress 
distribution along an axis of symmetry, perpendicular to the axis of 
rotation. The graphs are plotted from right to left. The radial and tan- 
gential stresses are maintained at the maximum for almost 80% of the 
radius and decrease only at the tip. 

These techniques permit a very accurate determination of all stresses 
throughout a homogeneous flywheel, and provide all the quantitative in- 
formation necessary to perform sensitivity tradeott studies. This allows 
a flywheel to be geometrically o;^imized for a given application in a pre- 
cise, quick, and economical fashion. The desirability of a constant stress 
and homogeneous material was assumed; however, in manufacturing thick 
forgings the metallurgical composition can vary considerably from Uie 
core to the surface as well as from the hub to the rim. The effects of the 
resulting stress pattern variations, developed within the material, must 
be taken into account and applied to the optimisation procedure presented. 
Such methods have been developed at Lockheed Missiles Ic Space Company 
in the form of additional computer programs that evaluate the effects of 
non-homogeneity of the flywheel matertal. 
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Table 1. Flywheel shape factors for various geometries 


Flywheel Geometry 

Sliape Factor 

Constant-stress disc (OD — ^ qp) 

1 . 00 

Modified constant -st res s disc (typical) 

0. 931 

Truncated conical disc 

0. 806 

Flat unpierced disc 

0. 606 

Thin rim (ID/OD— ^ 1.0) 

0. 500 

Shaped bar (OD— 

0. 500 

Rim with web (typical) 

0. 400 

Single filament (about transverse axis) 

0. 3 33 

Flat pierced disc 

0. 305 

From Ref. 1. ^ . 
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Table Z. Flywheel geometry and stress distribution 


rLYWHK^L PARAM^TFRS 


SPEED 

’MATERIAL OENStTY ------ 

P0IS50M3 PATIO 

YOOMTS MODJUJS 

TEMP. EXP. COErriCIEMT 
TOTAL TEMP. DlPTEWENCE 
OUTSIDE RADIUS -------- 


nffc> PAD/SEC 
♦PAS LBS/CU-IN 
.30 

30000. 0E^03 
T300.0E-09 
.0 DEG-F 
e 1.0000 INCHES 


FLYWHEEL DESIGN FACTORS 


MEAN TANGENTIAL STRESS 
MOMENT OF INERTIA ---- 

kinetic energy ------ 


WEIGHT . — — — 
ENERGY DENSITY 


SHAPE FACTOR - 
RADIAL GROWTH 


S 3 inSTPl PSt 

400. A53 LE-IN-SEC*e 

386417803 IN-LSS 
16.0066 KWH 
108P.41 LP 
3S9954«08 IN-LBS/LB 
11*3 WH/LB 
.9880 

•463714E-01 INCHES 


RADIUS 

THICKNESS 

TEMPERATURE 

TANGENTIAL 

RADIAL 

<IN) 

f 1N> 

(0KQ-F> 

STRESS 

STRESS 




(PSD 

<psn 

61.0000 

1*6658 

• 0 

86646 

0 

60*5000 

1*6858 

• 0 

76139 

13118 

PO.OOOO 

1*6658 

• 0 

77789 

88085 

14.5000 

1*6858 

.0 

83184 

38789 

19.0000 

1*6858 

.0 

88314 

51840 

18.5000 

1*6858 

• 0 

93188 

63569 

lA.oono 

1*6858 

• 0 

97734 

75787 

17.5000 

1 * 6658 

• 0 

101999 

87788 

17.0000 

1*6858 

.0 

105948 

99589 

18.8000 

1*3486 

.0 

107789 

104940 

16.0000 

1*4914 

.0 

108083 

105571 

IS. 5000 

1*6487 

.0 

108885 

106188 

15.0000 

1*8164 

• 0 

108519 

106881 

14.5000 

1*9885 

.0 

108729 

107058 

14.0000 

6*1748 

.0 

108917 

107446 

13*5000 

6*3716 

• 0 

109088 

107791 

11*0000 

6*5776 

.0 

109838 

108099 

16*3000 

6*7963 

.0 

109375 

108373 

16*0000 

3*0158 

• 0 

109499 

108819 

11*5000 

3*6471 

• 0 

109810 

108838 

11*0000 

3*4850 

*0 

109711 

109035 

io*50or 

3*7687 

• 0 

109806 

109611. 

10*0000 

3*9769 

• 0 

109885 

109369 

9 .5000 

4*6683 

*0 

109959 

109511 

9*0000 

4*4815 

• 0 

110087 

109639 

8*5000 

4*7350 

• 0 

110088 

109753 

8*0000 

4*9870 

.0 

110144 

109858 

T*S000 

5*6380 

• 0 

110194 

109948 

T.nooo 

5*4806 

• 0 

110839 

110031 

8*5000 

5*7178 

*0 

110879 

110104 

8.0000 

5*9470 

• 0 

110318 

110170 

5*5000 

6 1859 

• 0 

110349 

1 10889 

:«noon 

8.3768 

*0 

110377 

110881 

4*5000 

8*5860 

• 0 

110403 

110388 

4.0000 

8*7438 

• 0 

110485 

110388 

3.5000 

8*9047 

• 0 

110444 

110400 

3.0000 

7*0476 

• 0 

110461 

110489 

t*5000 

7*1701 

*0 

110474 

: 10454 

6.0000 

7*6766 

• 0 

110484 

110474 

I .5000 

7*3567 

*0 

t 104V1 

110490 

1*0000 

7*4107 

•0 

110493 

110505 

• 5000 

7*4457 

• 0 

110485 

110543 

• 0000 

7*4574 

*0 
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